The double-layer of penetrable ions: an alternative route to charge reversal 
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We investigate a double-layer of penetrable ions near a charged wall. We find a new mechanism 
for charge reversal that occurs in the weak-coupling regime and, accordingly, the system is suitable 
for the mean-field analysis. The penetrability is achieved by smearing-out the ionic charge inside a 
sphere, so there is no need to introduce non-electrostatic forces and the system in the low coupling 
limit can be described by a modified version of the Poisson-Boltzmann equation. The predictions 
of the theory are compared with the Monte Carlo simulations. 

PACS numbers: 



I. INTRODUCTION 

Overcharging is a surprising phenomenon in which con- 
terions adsorbed to a surface exceed the number of fixed 
surface charges [H. As a consequence, colons are drawn 
from the bulk toward the overcharged surface leading to 
a formation of a "triple-layer". Because the underlying 
mechanism of overcharging relies on electrostatic corre- 
lations 0, it has become synonymous with the strong- 
coupling limit. Indeed, a mean-field treatment which 
accurately captures the weak-coupling limit, cannot de- 
scribe overcharging. (However, the mean-field in com- 
bination with the excluded volume interactions can in- 
duce overcharging if the bulk volume fraction of an elec- 
trolyte is sufficiently large to generate a depletion force 
that pushes particles against a surface. This effect is seen 
in uncharged systems and persists for weakly charged sur- 
face charges 

In the strong-coupling limit overcharging is the result 
of increased structuring within the layer of counterions. 
The electrostatic correlations between the condensed 
counterions lead to formation of correlation "holes" 
within the layer of condensed ions which can attract ex- 
cess of counterions from the bulk. The value of the cou- 
pling constant T = Z^Xb / \gGi which is the ratio be- 
tween the Bjerrum and Gouy-Chapman length {Z is the 
valency), estimates the extent of correlation effects. In 
the limit F — ^ oo, the counterions are said to freeze into 
a 2d Wigner crystal At large, but finite F, the lo- 
cal structure of an ionic fluid remains Wigner- like 0-113 ■ 
The above mechanism is specific to Coulomb interactions 
that diverge as r — >• and, therefore, exhibit the ex- 
cluded volume effects d, [U. If, however, the divergence 
in the pair interaction is removed (the pair potential is 
bounded as r — > 0), particles can interpenetrate and the 
usual excluded volume interactions underlying the crys- 
tal formation are eliminated (at sufficiently high tem- 
perature and/or density). For some family of bounded 
potentials, particles can form stacks where two or more 
particles occupy the same position and act as an effective 
single particle. This is possible only if a pair potential is 
sufficiently fiat around r = One example is the 

penetrable sphere model where the pair potential is the 



step function Il3l4l6l . The stacking formations stabilize 
the liquid phase [ij, since doublets, triplets, etc. effec- 
tively decrease the number of particles. The presence of 
stacked formations is signaled as a positive peak in the 
pair correlation function at r — > Extrapolating 
these ideas to ions, which in addition to penetrable cores 
have long-range Coulomb interactions, we ask what influ- 
ence penetrability has on the structure of a double-layer. 
Can the restructuring invoked by the penetrability lead 
to overcharging in the weak-coupling limit? 

In the present work, the divergence in the Coulomb 
potential is removed by smearing-out the central charge 
of an ion over a finite region. The penetrating core, then, 
depends on the weight function used to smear out the 
charge. This procedure docs not require going beyond 
the framework of electrostatics and the weak-coupling 
limit can be described by the modified version of the 
Poisson-Boltzmann equation. 

Bulk properties and phase diagram of penetrable ions 
have been investigated in jF7H19j . The main feature is 
the formation of the Bjerrum pairs of two opposite ions 
that function as polarizable particles. The formation of 
these pairs leads to an insulator-conductor transition jT?! , 
which docs not exist in systems of hard-core ions in three 
dimensions [20j . Thus, penetrability dramatically affects 
the phase transition and the topology of the phase dia- 
gram. 

The model of penetrable particles is not only of the- 
oretical interest. Various macroparticles can exhibit in- 
terpenetration. Marquest and Witten [l^ suggested a 
penetrable sphere model for micelles. Polymer coils and 
dendrimers in good solvent can be represented by a Gaus- 
sian core model [2l|, . If in addition these macropar- 
ticles are charged, as is often the case for real systems 
[23| . then the model of penetrable ions can be of gen- 
uine physical relevance. Recently, ionic mic rog els have 
been modeled as uniformly charged spheres j24| . allow- 
ing interpenetration at short- and Coulomb interactions 
at long-separations. 
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II. THE POISSON-BOLTZMANN EQUATION 
FOR SMEARED-OUT CHARGES 



The charge qi of an ion i is smeared out over the finite 
region according to the weight function Wi(r — r^) such 
that / drwi(r — Ti) ^ qi. The charge density operator for 
N smeared-out ions is 

N 

Pc{r) - '■»)^»(^' - 

i—l 

In this work wc consider a symmetric electrolyte, qj^ ~ 
—q^ . The averaged charge density for this system is 



Pc(r) 



dv' 



p+ (r')w+ (r' - r) + p_ (r')w_ (r' - r) 



, (2) 



where pj^ and p_ denotes the number density of anions 
and cations, respectively. The Poisson equation is 



dv' 



P+{y')uj+{y' - r) + p_(r')c^-(r' - r) 



(3) 

where e is the background dielectric constant. To obtain 
a closed equation, we need an expression for p_ and p^ 
in terms of the mean electrostatic potential "0. For point 
charges this leads to 

p„(r) = c,e-^«°'^W, (4) 

where the subscript a is either + or — , and Cs is the 
bulk salt concentration. The number density depends 
locally on an electrostatic potential. However, if charge 
is smeared around the ion center at r, the entire dis- 
tribution uj{y' — r) interacts with the mean electrostatic 
potential, 



drV(r')wa(r'-r), 



and the number density becomes 

n =r p-/3/dr'V>(r')c.„(r'-r) 



(5) 



(6) 



We may now write down the mean-field equation for the 
electrostatic potential produced by smeared-out ions. 



eV'^ip = Cs dr'w+(r' - r)e 



-/3/dr"V>(r")L^+(r"-r') 



Cs J (ir'cj_(r' — r)e 



-/3 / dr"i/'(r")'^-(r"-r') 



(7) 



We refer to this modified Poisson-Boltzmann equation 
the Finite-Spread PB equation (FSPB). The FSPB equa- 
tion complements the already quite sizable set of modi- 
fied PB equations: the PB that incorporates the excluded 
volume interactions [25|, the dipolar interactions [26j . 



the nonlinear solvent contributions [27[, and the polar- 
izability of ions [2^, [2^ . The idea of finitely spread-out 
ions was considered in [s^l to study the self energy con- 
tributions beyond the mean-field. Mathematically, the 
FSPB equation has the structure of a convoluted equa- 
tion. Each ion is convoluted according to the weight func- 
tion which determines the composition of a single ion. 

The Eq. ([7]) can also be obtained by minimizing the 
grand potential. 



^^'^"^^ 2./ ./ ATier'-r 



kBTj2 / rfrp„(r)[logp„A3-l 
J drpa{r), 



(8) 



where A is the thermal de Broglie wavelength and p. = 
= /i_ is the chemical potential. The three contribu- 
tions are the electrostatic energy, entropy, and we allow 
the number density to fluctuate through the contact with 
a reservoir. The minimization = recovers Eq. © 
and the application of the Poisson equation yields the 
FSPB model. 



HI. VARIOUS DISTRIBUTIONS a;(r) 

The concept of charge smearing is not novel to this 
work, but it has been evoked many times in the past 
both as a physical representation and a mathematical 
construct. The best known example (and to our knowl- 
edge the earliest) of mathematical construct is the calcu- 
lation of Ewald summation to treat periodic charges [lH , 
today the most practiced method to account for contribu- 
tions due to periodic boundary conditions of a simulated 
system [s^], where by spreading the charge one achieves 
separation of interactions into the short- and long-range 
counterpart. Another instance of mathematical conve- 
nience is the Onsager smearing optimization to obtain 
the variational free energy of the strongly correlated one 
component plasma [s^. The idea of charges which at 
short separations exhibit soft repulsion and at long sepa- 
rations the Coulomb interaction appeared in [3Jj to rep- 
resent semi-classical hydrogen plasma at high tempera- 
ture and density, yet no smearing procedure was used 
to construct this potential. The actual smearing proce- 
dure to represent physical particles was used to model 
electrons set in cyclotron motion by a uniform magnetic 
field [H, H^. In soft matter, the smearing out procedure 
was used to represent microgels in (37| . and recently to 
represent charged polymers |17H19| . 

We consider a few simple, spherically symmetric distri- 
butions uj{r) and the interactions that result from these 
distributions. A spherically symmetric charge distribu- 
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tion generates the following electrostatic potential, 



ds s uj{s) + / dssu;{s) 



(9) 



To derive the above formula we applied the Gauss law to 
the Poisson equation, and afterwards we integrated the 
resulting electrostatic field (the integration by parts was 
evoked). If the distribution is uniform inside a spherical 
cavity. 



uj{r' — r) = 



3Q 



9{R-\r' -r\ 



(10) 



where 9 is the Heaviside function and / dr'uj{r' — r) = Q, 
then the potential inside a sphere is 



IV. NUMERICAL RESULTS FOR THE FSPB 
EQUATION 

We primarily focus on the distribution 

w±(|r - r,|) = ±ee{R - |r - r,\)/iy. (19) 

where i' = 47ri?'^/3 is the ionic volume. The FSPB equa- 
tion in reduced units reads 

.2,-1 / j„/ 



V^^ = K^i^-^ / dr'0{R - |r' - r|) sinh(0e), 



where 



e(r) = - / drV(r')^(i?-|r'-r|), (20) 



iP{r < R) 



Q_ 

Aire 



3R^ - 
2i?3 



(11) 



Outside the sphere the Coulomb potential is recovered. 
The interaction between two ions with this charge distri- 
bution is in 



U{r < 2R) 



91. 



192i?5 - SOR^r^ + 30R^r^ - 
160i?6 



(12) 

For non-overlapping separations, ions behave like point 
charges. Next, we consider the distribution 



^(r'-r) = ^<S(i?-|r'-r| 



(13) 



where a charge is smcarcd-out over a spherical shell. In- 
side a sphere, the electrostatic potential is constant. 



(14) 



The interaction between two of these distributions for 
overlapping separations is 



C/(,.<2R).|i 



4i?-r 
4i?2 



A sphere and a shell distributions interact via 



C/(,.<2fl).J- 



16i?4 



(15) 



(16) 



Finally, we consider a Gaussian distribution 



uj{r — r) 



(27r)3/2i?3 

which leads to the following pair interaction 

Q2 erf(r/2i?) 



U{r) 



47re 



(17) 



(18) 



(j> = Peip, — l/y/SirXsCs is the screening length, 
and Xb = iSe^/Ane is the Bjerrum length. The number 
density of each species is 



T<i>B 



(21) 



All the ions occupy the half-space x > 0. The charged 
hard wall at x = — i? determines boundary conditions, 



dx 



(22) 



where Cc is the surface charge. Note that the number 
density is limited to the region x > 0, but the charge 
density, because of the distribution ui{r), reaches all the 
way to X = —R, the location of the charged wall. Con- 
sequently, we assume the (hard-sphere)-(hard-wall) po- 
tential between an ion and the wall. Subsequent figures 
show data points for x > 0, the region available to ion 
centers. 

In the dilute limit (c^ — > 0), overcharging does not hap- 
pen, the solution of the PB equation for point ions in the 
weak-coupling limit yields an algebraic density profile. 



p{x) 



2t:Xb(tI 



(1 



(23) 



In the strong-coupling limit the electrostatic correlations 
modify the functional form of the distribution [38| . 



p{x) — 2TTXB<yl cxp( 



^^Gc)' 



(24) 



where Aqc — (SttAbCc)"^ is the Gouy-Chapman length. 
The barometric-like distribution in the transverse direc- 
tion is a consequence of high degree of ordering in the 
lateral direction reminiscent of the Wiener crystal 

Although the FSPB equation is purely mean-field, we 
find similar modification of the density profile, but this 
time the relevant parameter is R. The limit i? — > corre- 
sponds to the ionic distribution given by Eq. ((23)) . and the 
limit i? 00 to the one given by Eq. (P^ . see Fig. ([1]). 
In the limit i? — )■ 00, see Eq. (|T^ . we recover an ideal gas 
particles in a uniform gravitational field. 
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FIG. 3. The coion distribution at a charged wall correspond- 
FIG. 1 The counterion density profile in the dilute limit (no -^^ ^-^^ ^^^^^^ pjg q ^^-^^^ ^^^^^^ 

coions^The dotted lines demarcate the limiting behaviors in ^eutrahze the inverted charge. Note that the density profile 
Eq. m and Eq. The wall charge is positive and the exhibits oscillations, 

counterions have negative charge. As R increases, the profile 
takes the exponential shape. 



Fig. ([5]) shows the counterion profiles for an electrolyte 
solution. This introduces a length scale k~-^. Coun- 
terions from the bulk can now overcharge the surface. 
The counterion profile beyond x = 0.5nm dips below 
the bulk value, which indicates an overcharged surface. 
Fig. dSl) shows how the coion density rises above the 
bulk value, another signature of charge reversal. Finally, 
Fig. (jH) shows the non-monotonic electrostatic potential 
which goes to negative values and has a minimum around 
X = 0.6nm, at which point the electrostatic field vanishes 
and, farther on, it changes sign. 



— R=0 

— R=0.8nm 




\\ I I I I I i 1 1 T H T--T T ! 1 

0.5 1 1.5 

x[nm] 



FIG. 2. The counterion density profile near a charged wall. 
R — corresponds to point ions. The relevant lengths are: 
= 0.3nm A_B = 0.72nm, and Agc ~ 0.09nm. 




0- 



~' 1 2 3 

x[mn] 

FIG. 4. Electrostatic potential near a charged wall corre- 
sponding to the system in Fig. ((2]). The minimum corresponds 
to a point of vanishing electrostatic field. 



the thermal fiuctuations must exceed the interaction en- 
ergy of the overlapped particles, /3U{r = 0) ^ 1, or the 
counterion concentration at the wall must be sufficiently 
large to involve numerous overlaps, such that the effective 
2D density within the counterion layer is P2dR^ ^ 1- For 
the distribution in Eq. (|19p the first requirement translate 
into R/Xb ^ 1.2, and, assuming p2D ~ fc the second 
one becomes acR^ ^ 1. In Fig. ^ R/Xb ~ 1-1 and 
R^cTc ~ 1.6, and so both criteria are met. 

As a parenthetical note, we point out that the PB and 
the FSPB models satisfy the same contact value relation, 



Under what conditions penetrability of ions leads to 
charge reversal? For penetration to take place, either 



kBTpw = -Pid + TT' 



(25) 
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where Pid = ksTph is the ideal gas pressure in the bulk, 
Pb and is the total density in the bulk and at the 
contact with a wall, respectively. 

A. Tuning of the short-range interactions 

The sole constraint that a;(r'— r) needs to satisfy is that 
its integral recovers the charge of an ion. This leaves suf- 
ficient room to manipulate the short-range interactions. 
As an example we can take a mixed distribution 



w±(|r - r^l) = ±e 



{Z + l)0{R-\r-r^\)/i^ 



Z<5(i?-|r-r,|)/7 



(26) 

comprised of a charged shell and sphere. 7 = AttR^ is the 
surface area of a sphere. Z is the additional parameter 
and Z = corresponds to the distribution in Eq. ([T^. 
The number density is 



T [(2 + 1)08-205] 



where 



Mr) 



dr>(r')(5(i? 



(27) 



(28) 



Fig. ([5]) plots the pair potentials for different Z. The 
parameter Z controls the strength of the repulsion and, 
by the same token, the strength of the excluded volume 
interactions. By way of example. Fig. ^ shows the den- 




1 

r/R 



FIG. 5. Pair potential for two overlapping identical ions with 
charge distribution in Eq. (|26l) . Z — corresponds to the 
distribution in Eq. (fTOjl . 



sity profile for Z — and Z = 4, for R = 0.8nm. The ex- 
cluded volume contributions for Z — A expel counterions 



from the first layer, which is the opposite of overcharging 
seen for Z = 0. 

With this example we try to demonstrate possible ex- 
tensions of the model based on charge spreading. For 
example, the repulsion controlled by the parameter Z 
could represent the interaction between polymers on ac- 
count of the self-avoiding walk of polymer chains. 




z=o 

Z=4 



x[nm] 



1.5 



FIG. 6. The counterion density profile near a charged wall for 
the point ions and for the the ions with the distribution in Eq. 
(|26|) for R — O.Snm and different Z. The system parameters 
are those in Fig. ([2]). 



V. COMPARISON WITH SIMULATION 

To test the mean-field approximation implicit in the 
FSPB equation, we carry out Monte Carlo simulations. 
The size of the simulation box is Ly — ~ Sum and 
Lx = 6nm. The periodic boundary conditions arc in the 
lateral (y, z)-directions. The bounding plates at x' = 
and X = Lx have opposite charge. The box encloses 
N = 600 ions, N+ = = 300. To check for the finite 
size effects, we double the simulation size for some pa- 
rameters. The standard Ewald summation is applied for 
contributions of periodic images [32| . If particles over- 
lap, we supplement the interaction energy with the term: 
l3AUtot{r < 2R) = l3U{r,j) - (A physical picture is 
slightly modified when doing simulations. In simulations 
we tend to think of particles as point charges which at 
separations r < 2R switch to the non-Coulomb interac- 
tions U{r). Within the mean-field theory based on the 
Poisson equation, we tend to think of an ion as a charged 
ball. The two conceptions arc, however, identical. ) 

Fig. ([7]) compares counterion profiles obtained from the 
simulation and the FSPB equation. For R = O.Snm the 
agreement is virtually exact. For the smaller radius R = 
O.lnm we see the two results deviate: in the simulation, 
the correlations generate a more concentrated layer of 
counterions (but no overcharging). In Fig. ([S]) we plot 
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the coion density profile for R = 0.8nm to further confirn] 
the accuracy of the mean-field in this regime. 




x[nm] 



FIG. 7. Counterion density profiles. The system parameters 
are: Agc ~ 0.09nm and Xb = 0.72nm. The symbols designate 
the simulation data points and the lines are the numerical 
results for the FSPB equation. 




x[nm] 



FIG. 8. The coion density profile for R = 0.8nm obtained 
from the simulation (symbols) and the FSPB equation (solid 
line). The remaining parameters are as in Fig. ((T)). 

What causes charge reversal for smeared-out ions? Do 
counterions merely penetrate or the pair potential is sufR- 
ciently flat in the overlap region and counterions collapse 
into stacked formations, so that the pair correlation func- 
tion exhibits positive peak as r 0? In Fig. ([9|) we show 
the configuration snapshots for counterions adsorbed on a 
charged surface for different values of R. The difference in 
structure is perceptible. Ionic penetration favors smaller 
separations between counterions. This gives the impres- 
sion of string-like formations. As mentioned in the intro- 
duction, not all penetrable potentials can lead to stacked 
configurations, where two or more particles occupy "the 
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FIG. 9. Configuration snapshot of counterions adsorbed on 
the charged wall at a; < 0.35nm. The diameter of particles is 
a = 0.5nm and is selected arbitrarily for visualization. The 
configuration on left is for R/Xb ~ 0.14 and the one on the 
right is for R/Xb = 1.1- The 2d densities are p2d = 2.34nm~^ 
and p2d = 2.68nm~^, respectively. For comparison, the sur- 
face charge density is ctc — 2.50nm^^. The system parameters 
are as in Fig. Q. 



same" space. For example, the Gaussian core model pre- 
cludes [3^ , while the penetrable sphere model favors 
the stacked configurations [13, [3- Stable stacked ag- 
gregates require that a pair potential be sufficiently flat. 
A more rigorous test involves the Fourier transform of 
the pair interaction. If U{k) oscillates, and, therefore, 
involves negative values, then the stacking takes place 
under certain conditions [l2[. On the other hand, if U{k) 
is non-negative, stacking docs not occur under any con- 
ditions. This criterion is arrived at by considering the 
Ornstcin-Zernike relation which, in the Fourier space, is 



hik) 



c(fc) 



1 - pc{k) ' 



(29) 



where c(r) and h{r) is the direct and pair correlation 
function, respectively. In the mean-field approxima- 
tion cmf(^) ~ -'l3U{r) [the exact definition is c(r) = 
— sp(^JIs'p^j.,^ ■ a-nd in the mean- field F^^. is the first term in 
Eq. dH])], so that 



pUjk) 
l + pl3U{k)' 



(30) 



If we take h{r = 0) > to be the signature of stacking, 
and 



(31) 



then it becomes immediately clear that a non-negative 
U{k) always yields /imf(O) < 0, and no stacking oc- 
curs. Only an oscillating U{k) can yield /imf > 0- In 
the present model particles are smeared-out charges. 



pU{\v~r'\) = Xb / dr" / r' 



uj{v - r")w(r' - r'") 



(32) 
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and in the Fourier space 



/3f/(fc) 



47rABa;^(fc) 
fc2 ■ 



This yields 
/imf(O) 



1 

2^ 



dk 



fc2 + 47rAB/9w2(fc) ' 



(33) 



(34) 



We see that h{r = 0) < regardless of the distribution 
U!{r). We conclude that soft interactions generated by 
charge spreading cannot lead to stacked configurations. 

In Fig. (jlOp we show the lateral correlation function 
for the adsorbed coutnerions, h\\{r). Penetration reduces 
the degree of electrostatic correlations between the ions, 
however, h\\{r) always remains negative and decreases 
monotonically all the way to r = 0. The small difference 
between the shell and the sphere distributions indicates 
that the exact shape of the pair potential is unimpor- 
tant. In Fig. prj) we show the probability distribution 



•— • sphere 
--shell 
Gaussian 




0.5 

r[nm] 

FIG. 11. The nearest neighbor separation distribution for 
different cj(r) functions. R — 0.8nm and the same parameters 
as in Fig. 0. 



-0.5 



• R=0.1nm (sphere) 
R=0.8nm (sphere) 
=.R=0.8nm (shell) 



r[nm] 



0.5 



FIG. 10. Correlation function for counterions in the lateral 
plane adjacent to the wall. The layer thickness is 0.35nm. 
The remaining parameters are the same as in Fig. ((T]). The 
dashed lines guide the eye. 

of the nearest neighbor separation for various Lu{r) func- 
tions. Differences are rather small and only quantitative. 
The Gaussian uj{r) shows preference for smaller interionic 
separation. 

Finally, we explore the validity of the mean-field for 
multivalent symmetric ions, Z > 1. We suppose that the 
mean-field should deteriorate quickly for Z > 1 since the 
coupling constant T ^ Z^ . Fig. (|T2|) shows the counterion 
density profiles for Z = 2,3. For Z = 2 the mean-field 
is virtually exact and for Z ~ i there are small devia- 
tions. This surprising agreement can be explained with 
the following. Increasing Z naturally reduces the num- 
ber of counterions required to neutralize the wall. On the 
other hand, a larger Z generates stronger overcharging so 
the density drop is partially compensated. The number 



of overlapping configurations is still large and the mean- 
field retains its validity. To confirm this conjecture, we 
check the 2D density of the counterions adsorbed on the 
charged wall. For Z = 3, we find P2D ~ 1-2 > 1. The 
overlapping configurations are, therefore, still significant. 




x[nm] 



FIG. 12. Counterion distribution function for symmetric so- 
lution with ion charge ±eZ. The remaining parameters are 
the same as in Fig. ([7|. The circles are the simulation data 
points and the lines are obtained from the FSPB. 



VI. DISCUSSION 

For overcharging to occur, there must exist a mecha- 
nism that either reduces the free energy of counterions 
condensed on a charged surface (it becomes favorable for 
a counterion to leave its bulk environment) or increase 
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the free energy of a bulk (it becomes less favorable for 
ions to stay in the bulk). It is the exchange between the 
two environments, the bulk and the double-layer, that 
leads to overcharging. We consider the environment of a 
charged surface with its adsorbed counterions and con- 
sider the conditions that lead to lowering of its energy. 
The surface charge of a plate is ■ We do not violate the 
condition of neutrality, thus the 2D density of adsorbed 
counterions is adZ. We focus on the energy of a single 
counterion. For sake of simplicity, we assume counterions 
and the charged wall occupy the same plane. An attrac- 
tion to the wall yields a negative energy contribution, 



wall 



lim 



lim 



27r 



dr r- 



— 2'KZXB(yci'c 



(35) 



which is countered by the repulsive interactions with 
other counterions, 



lim 



27r 



dr rg(r) 



ZXbct. 



(36) 



where g{r) is the pair distribution function. For absence 
of correlations g{r) = 1, and the attractive contribution 
is completely canceled out. 



lim 



iTlZXljCf^ 



(37) 



But the cancellation will not be exact if correlations are 
present. Separating correlations from the pair distribu- 
tion function, g{r) = /i(r), the uncanceled correlation 
energy is 



1nXB<yc 



dr h{r), 



(38) 



The relation between the correlation energy and the cor- 
relation function is quite clear. An ion will generate 
a negative correlation hole in its neighborhood, which 
will lead to negative energy contributions. We assume 
the following simple correlation hole, h{r > ^) = and 
h{r < Q) = —1, that is, by fixing a position of an ion, 
we introduce a circular hole in the density profile with 
the correlation length C,. The conservation of mass con- 
dition requires that an area of the hole is related to the 
density via adZ = l/(7rC^), therefore, the correlation 
length is C = \/ Zj (ttCTc)- This correlation function is the 
low temperature limit of the correlation hole obtaine d by 
construction from the linear mean- field treatment 41 1. 
Inserting this h(r) into Eq. (j38p we can approximate the 
correlation energy in the strong-coupling limit. 



(39) 



Recalling the definition of the coupling constant F 
27rZ^A^ (Tc, we get 



'2Y. 



(40) 



But for penetrable ions correlations are not required 
for reducing the electrostatic energy of counterions - the 
fact demonstrated by the validity of the mean-field the- 
ory. Consequently, we set h[r) = 0. The reduction in 
electrostatic energy comes from a different source, from 
the fact that at overlapping separations, r < 2R, the elec- 
trostatic interactions are reduced on account of smearing 
procedure, which leads to the energy gain 



/3 Overlap « 27r / dr Ocr 
Jo 

= —CttZXbcTcR, 



/3U{r) 



ZXf 



(41) 



where the constant C depends on the pair potential U{r), 
which, in turn, is determined by the distribution func- 
tion Lo{r). For the distributions considered in this work: 



sphere 



36/35, a 



shell 



4/3, and Cg 



2.06. Com- 



paring Eq. ([M]) with Eq. (|1T|) we see the different depen- 
dence of each mechanism on different parameters. The 
stabilization based on penetration has stronger depen- 
dence on the surface charge and the Bjerrum constant, 
on the other hand, its dependence on valency is weaker. 
In the strong-coupling limit, penetrable ions will exhibit 
correlated motions. If the correlation length is larger 
than the diameter of an ion, ^J'zjijxad) > 2i?, penetra- 
tion may be neglected and the forme r mechan ism comes 
to the fore. On the other hand, if Z/ (ttcTc) < 2i?, we 
expect the two contributions to mix. 

To recap, both mechanisms depend on eliminating 
the energy contributions coming from short separations 
between ions. For point ions in the strong-coupling 
limit, configurations with short separations are elimi- 
nated through correlated motion. The price is sacrifice 
in entropy, despite this, the total free energy is lowered. 
For penetrable ions the problem of high energy contribu- 
tions at short separations does not exist to begin with. 
Due to smearing out procedure of an ion charge and the 
removal of the divergence from pair interactions, these 
contributions are taken out of the picture. Consequently, 
there is no entropy price to be paid, as all separations 
are explored " equally" and the mechanism is valid in the 
weak-coupling limit. 



VII. CONCLUSION 

The present work studies the structure of a double- 
layer composed of ions whose central charge is smeared 
over a finite region in accordance with a weight function 
u>{r). The smearing-out procedure removes the diver- 
gence as r — > from the pair interaction, allowing for 
interpenetration between the ions. The conditions under 
which penetration is favored are large temperature and 
high density. This regime is suitable for the mean-field 
treatment. Accordingly, we derive a modified Poisson- 
Boltzmann equation for spread-out charges (the FSPB 
equation). The FSPB equation predicts that for suffi- 
ciently large spreading radius R, overcharging takes place 
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- the MC simulations verify this prediction. This sug- 
gests an ahernative mechanism for charge reversal that is 
not related to correlations and the strong-coupling limit. 

Using simulations and the mean-field approximation, 
we can exclude the ionic stacking as the underlying mech- 
anism of overcharging. In fact, any soft repulsion ob- 
tained by charge spreading cannot lead to stacked con- 
figurations. 

As a final consideration, we address the physical rel- 
evance of the smeared-out ion model. The spreading- 
out of the charge may capture the interactions between 
charges distributed along the polymer chains, but a more 
realistic representation would involve a non-electrostatic 
component produced by the self-avoidance of the poly- 
mer chains. Within the electrostatics framework, we have 
suggested a plausible distribution ui{r) composed of a 



charged shell and sphere, which generates an additional 
repulsion inside a penetrable core. In the end, however, a 
physically accurate pair potential requires corroboration 
with experimental analysis. 
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